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We study the linear sigma model with O(N) symmetry at nonzero temperature in the framework 
of the Cornwall-Jackiw-Tomboulis formalism. Extending the set of two-particle irreducible diagrams 
by adding sunset diagrams to the usual Hartree-Fock contributions, we derive a new approximation 
scheme which extends the standard Hartree-Fock approximation by the inclusion of nonzero decay 
widths. In this approximation, in-medium modifications of the meson masses as well as decay and 
collisional broadening effects are self-consistently taken into account. As compared to the standard 
Hartree-Fock approximation, the temperature for the chiral symmetry restoring phase transition is 
lowered by ~ 20 — 25%. For large temperatures, the spectral densities of the cr-meson and the pion 
become degenerate, as required for the chirally symmetric phase. 


I. INTRODUCTION 

For zero quark masses, the QCD Lagrangian exhibits a global chiral U(Nf) r x U{Nf)z or, equiv¬ 
alently, U{Nf)v x U(Nf)A symmetry, where Nf is the number of quark flavors, and V = r +1, A = 
r — i. In the vacuum, the U(1)a anomaly [lj breaks this symmetry explicitly to U(Nf)y x SU(Nf)A- 
Furthermore, a nonvanishing chiral condensate (qq) ~ (300 MeV) 3 breaks this symmetry sponta¬ 
neously to U(Nf)v !, leading to Nj — 1 pseudoscalar Goldstone bosons which constitute the 
effective degrees of freedom at low energies. In nature, chiral symmetry is also explicitly broken by 
nonvanishing quark masses, which accounts for the physical masses of the pions, kaons, etc. 

At temperatures of the order of ~ ( qq ) 1//3 , the thermal excitation energy is large enough to expect 
the restoration of chiral symmetry. At such energy scales, the QCD coupling constant is still large, 
rendering perturbative calculations unreliable. Thus, one has to resort to nonperturbative methods 
to study chiral symmetry restoration. A first-principle approach is lattice QCD 53]. Lattice QCD 
calculations have determined the temperature T c for chiral symmetry restoration to be of order 
150 MeV at zero quark chemical potential 0. These calculations, however, face several technical 
problems. 

The first is that they become numerically very difficult for physically realistic, i.e., small, values 
of the up- and down-quark masses. Although progress in this direction has been made J5|, most 
studies use unphysically large values. Another problem is that, at nonzero chemical potential, lattice 
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QCD calculations are hampered by the fermion sign problem and become increasingly unreliable for 
chemical potentials larger than (a factor 7r times) the temperature Jb). 

An alternative nonperturbative approach to study chiral symmetry restoration is via chiral effec¬ 
tive theories. These theories have the same global U(Nf) r x U(Nf)e symmetry as QCD but, since 
quark and gluons are integrated out, do not possess the local SU(3) C color symmetry of QCD. The 
effective low-energy degrees of freedom are the (pseudo-) Goldstone bosons of the QCD vacuum, i.e., 
the pseudoscalar mesons. However, in the chirally symmetric phase these particles become degener¬ 
ate with their chiral partners, the scalar mesons. Therefore, an appropriate effective theory to study 
chiral symmetry restoration in QCD is the linear sigma model yj [8] which treats both scalar and 
pseudoscalar degrees of freedom on the same footing. The advantage of chiral effective theories over 
lattice QCD calculations is that their numerical treatment (within some many-body approximation 
scheme) is comparatively simple and that there is no problem to consider arbitrary quark chemical 
potential. 

At nonzero temperature T, ordinary perturbation theory in terms of the coupling constant g 
breaks down in the sense that one can no longer order different contributions according to powers 
of g 0]. This is because the new energy scale introduced by the temperature can conspire with the 
typical momentum scale p of a process so that gT/p is no longer of order g , but can be of order 1 
(lCL llH. Consequently, all terms of order gT/p have to be taken into account which requires the 
resummation of certain classes of diagrams. 

A convenient resummation method is provided by the extension of the Cornwall-Jackiw-Tomboulis 
(CJT) formalism 0 to nonzero temperatures and chemical potentials. The CJT formalism is 
equivalent to the ^-functional approach of Luttinger and Ward |l3j and Baym 0- It generalizes 
the concept of the effective action T[</>] for the expectation value <fi of the one-point function in the 
presence of external sources to that for the effective action 

T[</>, G] = S[fi + \ Tr In G" 1 + ± Tc{G~ l G - 1) + T 2 0, G] (1) 

for (j> and the expectation value G of the two-point function in the presence of external sources. 
Here, 5[0] is the tree-level action, G _1 the inverse tree-level two-point function, and T 2 [^, G] the 
sum of all two-particle irreducible (2PI) vacuum diagrams with internal lines given by G. (For an 
extension of this approach to three- and more-point functions, see [H,0ill3!la-) The stationary 
points of this functional, 

sma\ sr&G] 

^=ip,e=g ^G 

provide self-consistent equations for the expectation values of the one- and two-point functions cf> 
and G in the absence of external sources, denoted as (p and Q, respectively. The latter is the 
Dyson-Schwinger equation, 



g~ l = G ~ 1 + n, 


(3) 


where 


n = -2 


6TS,G] 


SG 


p=cp,G=g 


(4) 


is the self-energy. As long as T 2 contains all 2PI diagrams, the CJT effective action is exact. Of 
course, it is practically impossible to compute all 2PI diagrams, and one has to truncate T 2 , for 
instance at some order in the number of loops. The advantage of the CJT formalism is that any 
truncation of T 2 yields a many-body approximation scheme which preserves the symmetries of the 
tree-level action, if the expectation values cf> and G transform as first- and second-rank tensors of 
these symmetries 0]. The solution of Eqs. m is thermodynamically consistent and conserves the 
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Noether currents. However, Ward-Takahashi identities for higher-order vertex functions may be 
violated because of the truncation of T 2 [lj|. 

The most popular among these many-body approximation schemes is the Hartree-Fock approxi¬ 
mation. In this case, To solely contains diagrams of so-called “double-bubble” topology, cf. Fig. |T| 
a-c. Note that, in Refs. |2^[2i],[22j,|23j], this scheme has been termed “Hartree” approximation. How¬ 
ever, since the exchange contributions with respect to internal indices are in fact included, it is more 
appropriate to call it “Hartree-Fock” approximation. Neglecting the exchange contributions leads to 
the actual Hartree approximation which, in the case of the O(N) model, has also been termed “large- 
N ” approximation [2C|,[2l]|. In the Hartree-Fock approximation, the self-energies, which according 
to Eq. are obtained by cutting lines in the diagrams for F 2 , only consist of “tadpole” diagrams, 
cf. Figs.0 b,c, ancl[3]a,c. For the chiral effective theories, such as the U(Nf) r x U(Nf) f and Q(N ) 
models, this approximation scheme has been studied in great detail |20t I21L 12211231 124125L 126}. I27f. 
The Hartree-Fock approximation is a very simple approximation scheme, since tadpole self-energies 
do not have an imaginary part and, consequently, all particles are stable quasiparticles. Moreover, 
since tadpoles are independent of energy and momentum, the Dyson-Schwinger equations for the 
propagators reduce to fix-point equations for the in-medium masses. 

There are, however, several problems with the Hartree-Fock approximation related to the trunca¬ 
tion of F 2 . For instance, in the case of the O(N) model, it does not correctly reproduce the order of 
the chiral symmetry restoring phase transition and it violates Goldstone’s theorem in the sense that 
the Goldstone bosons do not remain massless for nonzero temperatures 0 < T < T c . Several ways 
to remedy this shortcoming have been suggested. The simplest one is to neglect subleadin g co ntri- 
butions in 1 /TV, leading to the Hartree (or large- N, see discussion above) approximation (20l l21l. 
Here, the Goldstone bosons remain massless for T < T c and the transition is of second order, as 
expected from universality class arguments 2Sj. Another possibility to restore Goldstone’s theorem 
is via so-called “external” propagators 0, I29L (30i . [dTl |32]. These objects are obtained from the 
Hartree-Fock propagators by additionally resumming all diagrams pertaining to the Random Phase 
Approximation (RPA), with internal lines given by the Hartree-Fock propagators. Another way to 
restore the second-order nature of the chiral phase transition is to include so-called “two-particle 
point-irreducible” (2PPI) contributions to T 2 1331133, 35], 

In this work, we are not concerned with these formal shortcomings of the Hartree-Fock approxima¬ 
tion: we focus exclusively on the case realized in nature where chiral symmetry is already explicitly 
broken by (small) nonzero quark masses, such that the pion is no longer a Goldstone boson and 
assumes a mass M n ~ 139.5 MeV in vacuum (36|. Our goal in this work is to include the nonzero 
decay width of the particles in a self-consistent many-body approximation scheme. To this end, one 
has to go beyond the Hartree-Fock approximation and add other diagrams to T 2 which, upon cutting 
lines according to Eq. o, yield self-energies with nonzero imaginary part. The most simple of such 
diagrams, and the ones considered in the following, are those of the so-called “sunset” topology, cf. 
Fig. n]d,e, leading to the self-energy diagrams shown in Figs. 0d,e and EH d. 

The self-energies arising from the sunset diagrams in F 2 depend on the energy and the momentum 
of the incoming particle. Therefore, the Dyson-Schwinger equations no longer reduce to fix-point 
equations for the in-medium masses. Since the self-energies now have a nonzero imaginary part, 
implying a finite lifetime of the corresponding particles, the spectral densities are no longer delta- 
functions. It is therefore convenient to rewrite the Dyson-Schwinger equations for the propagators 
into a set of self-consistent integral equations for the spectral densities which has to be solved 
numerically. 

The paper is organized as follows. In Sec. DU we discuss the application of the CJT formalism to 
the O(N) model. In Sec. IHIl we discuss our new approximation scheme which extends the standard 
Hartree-Fock approximation by including nonzero decay widths. Numerical results for the in-medium 
masses and decay widths of cr-mesons and pions, as well as their spectral densities are discussed in 
Sec. II VI Section Ivl concludes this paper with a summary of our results. Technical details are deferred 
to the appendices. 

We denote 4-vectors by capital letters, X M = (xo,x), with x being a 3-vector of modulus x = |x|. 
We use the imaginary-time formalism to compute quantities at nonzero temperature. Integrals over 
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4-momentum A' M = (fco,k) are denoted as 


r 00 _ r u 

/ f{K) = T^2 ?2^V3 k) , (5) 

JK n=-oo J ' ’ 

where T is the temperature and uj n = 2 tt nT, n = 0, ±1, ±2,... are the bosonic Matsubara frequen¬ 
cies. Our units are h = c = ks = 1. The metric tensor is = diag(+, —, —, —). 

II. THE O(N) MODEL IN THE CJT FORMALISM 

In this section we discuss the application of the CJT formalism to the 0(N ) model. In Sec. IIVI 
we shall present the numerical results for the case N = 4. The Lagrangian of the O(N) linear sigma 
model is given by 

£(0) = d^cf) -U{4>) (6) 

where </> = (</>i,7r), with the first component 4> 1 corresponding to the scalar a -meson and the other 
components 7 r = (</> 2 ,..., 4>n) corresponding to the pseudoscalar pions. The function U(4>) is the 
tree-level potential, 


U{4>) = 4>+ 4>) 2 -Hfa . (7) 

The constant /i 2 is the bare mass and A > 0 is the four-point coupling constant. For /i 2 < 0 the 
O(N) symmetry is spontaneously broken to 0(N — 1), leading to N — 1 Goldstone bosons, the pions. 
The parameter H breaks the symmetry explicitly, giving a mass to the pion. The determination of 
the parameters fjr , A, and H as functions of the vacuum masses of cr-meson and pion, as well as the 
pion decay constant, is discussed in detail in Appendix iBl 

We assume translational invariance so that we may consider the effective potential V instead of 
the effective action 0 For translationally invariant systems, these two quantities are related via 

T[a,it,S,P} = -^-V{a,7r,S,P), ( 8 ) 

where ^3 is the 3-volume of the system, and a. iv, S, P are the expectation values of the one- and 
two-point functions for the scalar and pseudoscalar fields in the presence of external sources |l 2 j . 
We are interested in the case where these sources are zero, i.e., in the stationary points of F or V, 
cf. Eq. 0. Because the vacuum of QCD has even parity, the expectation values of the pseudoscalar 
fields are zero, 7 r = 0, and we shall simply omit the dependence of V on n in the following. 

Then, the effective potential for the O(N) model in the CJT formalism reads [ 2 H |23| 

V[a,S,P] = U(a) + i f [\nS- 1 (Q) + S~\Q- 1 a)S(Q)-l} 

J Q 

+ I [\nP- 1 (Q) + p-\Q;a)P(Q)-l]+V 2 {a 1 S,P} , (9) 

1 Jq 

where U(a) is the tree-level potential 0, evaluated at 0 = (a, 0,... ,0). The quantities S ' -1 and 
P -1 are the inverse tree-level propagators for scalar and pseudoscalar mesons, 

S~\K-a) = -K 2 +m 2 a (a ), 
p-\K;a) = —K 2 + m 2 (ct) , 


(10a) 

(10b) 
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where the tree-level masses are 

ml{a) = + (11a) 

+ 77^ ■ ( llb ) 

The functional Vi in Eq. © is related to P2 in Eq. © via T 2 = —(Slz/T)V 2 . The Hartree-Fock 
approximation is defined by restricting the sum of all 2PI diagrams, V2, to the three double-bubble 
diagrams shown in Figs. |T| a, b, and c, 


v 2 a [P} 

V 2 b [S} 
V 2 C [S,P} 


(N + 1)(N — 1) 


N 


1 2 


P(Q) 


4 


Uq 


S(Q ) 


2 (N — 1) 




’Q 


S{Q) / P(L) . 


(12a) 

(12b) 

(12c) 


As explained in the introduction, in order to include the nonzero decay width of the particles we 
have to go beyond the Hartree-Fock approximation by additionally including the sunset diagrams 
of Figs.^d and e. These diagrams have an explicit dependence on <7, 


V 2 d [a,S,P] = ±2 (tf-l)(^) J l J s{L)P{Q)P(L + Q), (13a) 

S] = \ 3! (jf} J S(L)S(Q)S(L + Q) . (13b) 

The complete expression for V2 is obtained by the sum of the contributions (I12al) (I13bl) . 

v 2 = v 2 a + + V 2 + V 2 d + V 2 . (14) 


The expectation values of the one- and two-point functions in the absence of external sources, a and 
S , V, are determined from the stationary points of V, 


5V 

5a 


= 0 


a=a,S=S,P=V 


SV 

Js 


= o, 


S=a ,S=S,P=V 


SV 

Jp 


= o, 


a=a,S=S,P=V 


leading to an equation for the scalar condensate a , 

4A 4A r 

H = » 2 a + -a 3 + -aj' [3S(Q) + (N-1)V(Q)} + 

where 


SV* 

5a 


(15) 


(16a) 


(T—0,S—S,P—P 


5V 2 
5a 


a=a,S=S,P=V 


2(N — 1) / / S(L)V(Q)P(L + Q) 
Jl Jq 


+ 3! / / S(L)S(Q)S(L + Q) 
Jl Jq 


(16b) 


and to the Dyson-Schwinger equations for the scalar and pseudoscalar propagators, 

S~\K-a) = S~ 1 {K-,a)+'£[K-a], 

V~\K-a) = P- 1 {K-a)+H[K-a]. 


(16c) 

(16d) 
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Here we introduced the self-energies of the scalar and pseudoscalar fields, 


with 


and 


E = E a + E b + E c + E d + E e , 

n = n a + n & + n c + n d + n e 


E a 

E b 

E c 

E d [A>] 

S e [K;cr} 

n a 

n 6 

n c 

n d [K-a] 

n e 


= 2 


6Vn a 


= 2 


ss 

syl 

ss 


= o, 


a=a,S=S,P=V 


„ SV 2 C 
2 — 2 - 
6S 

„ 8Y 2 d 
2 — i- 
SS 


a=(T,S—S,P—V 


a=a,S=S,P=V 


a=a,S—S,P—V 


= ¥ 3 / 0 ^»- 

= i^)‘ 2{N - r 'S a p[K -^ Q '’ 


= 2 


SVo e 


6S 


<r=cr 1 S=S ,P—V 

2 SV 2 a 


>Q 

(^) 3 ' 3! f s ( k -qw®> 


N- 

1 

SP 

2 


sv 2 b 

N - 

1 

SP 

2 


SV 2 C 

N - 

1 

SP 

2 


sv 2 d 

N- 

1 

SP 

2 


sv 2 e 


a—cr,S—S,P—P 


= V< W +D 


= 0 


a=a,S=S,P=V 


= S / «S(Q) 


a=a,S=S,P=P 


N 


a =t r,S=5,P=P 


JV-1 


/ 4Act 

= V^/ 
= o. 


4 [ V(K - Q)S(Q) , 
Jo 


S =C r,5=5,P=P 


(17a) 

(17b) 

(18a) 

(18b) 

(18c) 

(18d) 

(18e) 

(19a) 

(19b) 

(19c) 

(19d) 

(19e) 


The calculation of the self-energy contributions Eqs. m and m corresponds to opening internal 
lines in the diagrams of Fig.|U Employing this procedure to the double-bubble diagrams leads to the 
tadpole contributions for the self-energies of the er-meson (see Figs. Ob, c) and the pion (see Figs.[J| 
a, c). This defines the Hartree-Fock approximation. Additionally, the sunset diagrams of Figs, ^d, 
e lead to the one-loop contributions shown in Figs. Qd, e and Eld. The latter contributions depend 
on the energy and the momentum of the particles and lead to nonvanishing imaginary parts for the 
self-energies. The explicit calculation of the self-energies m and m is discussed in Appendix 1X1 


III. IMPROVING HARTREE-FOCK BY INCLUDING NONZERO DECAY WIDTHS 

On the basis of the general framework presented in the preceding section we now discuss in detail 
how the Hartree-Fock approximation can be improved by including the nonzero decay width of the 
particles. 

As explained in the introduction, we rewrite the set of Dyson-Schwinger equations lllficl and (fTTTdl) 
in terms of a set of integral equations for the spectral densities. In general, the spectral densities for 
the cr-meson and the pion are defined as 

Pa(u), k) = 21m S(u + it], k; a) , 

Pn(u), k) = 21m V(ui + ip, k;cr) . 


(20a) 

(20b) 
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If the imaginary parts of the self-energies are zero, as in the Hartree-Fock approximation, all particles 
are stable quasiparticles, i.e., their spectral densities become delta-functions with support on the 
quasiparticle mass shell, 

p<T(w,k) = 27rZ CT [w CT (k),k] {<5[w - w CT (k)] - <5[w +o; cr (k)]} , (21a) 

Ptt(w, k) = 27rZ 7r [u; 7r (k),k] -w x (k)] - 6[u; + av(k)]} , (21b) 

where cJ a ,n( k) are the quasiparticle energies for a -mesons and pions, defined by the positive solutions 
of 


k) = k 2 +ml(a) + Re E[w CT (k), k; cr] , 
w^(k) = k 2 +ml(cr) + Ren[w^(k), k; a] , 

Z a [u a { k),k] 

Z w [u n ( k),k] 


dS-\K;a) 


dk 0 

dT~ l (K-a) 


dk 0 


ko-Ua- (k) 
-1 

kr\=UJ^( kl 


( 22 a) 

( 22 b) 


(23a) 

(23b) 


are the wave-function renormalization factors on the quasi-particle mass shell. Since the real parts 
of the self-energies are even functions of energy, these factors are also even in Wcr i 7 r (k). In the 
Hartree-Fock approximation, the (real parts of the) self-energies do not depend on such that 


where 


LOcr(k) = yjk 2 + M 2 (cr) , w^(fc) = y/k 2 + M 2 (a) , 


M 2 {a) = ml(cr) + Re£(er) , M 2 (cr) = ml (a) + Re n(cr) , 


are the effective er-meson and pion masses, and 

1 


Z a [uj a {k)] = 


2 u]<j(k) 


Zir[uj n {k)} = 


2w 7r (fc) 


(24) 

(25) 

(26) 


For nonvanishing imaginary parts of the self-energies the spectral densities assume the following 
form: 


Ar(w,k) 

Pi r(w,k) 


2 Im £(w, k; a) 

[u 2 — k 2 — to 2 (cr) — Re£(w, k; a)] 2 + [ImS(u, k; a)] 2 
2 Imn(u, k; cr) 

[w 2 — k 2 — ml (a) — Ren(w, k; cr )] 2 + [Imn(u, k; <r )] 2 


(27a) 

(27b) 


The general calculation of the imaginary parts of the self-energies is discussed in the second part 
of Appendix^ It turns out that they do not depend on the direction of 3-momentum k, but only 
on the modulus k. The tadpole diagrams of Figs. 0 b, c (for the er-meson) and Figs. 0 a, c (for the 














pion) do not have an imaginary part. Thus [see Appendix 1X1 Eq. IA19I ] . 


ImS(w, k;a) = Im E d (w, k; a) + Im E e (w, k; cr) 

{A\rr \ ^ 1 1 r°° 

= 2 (2tt ) 3 k J dw i dw 2[l + /(wi) + /(w 2 )]<5(w-wi-W2) 

POO 

x / dq 1 q 1 dq 2 q 2 0 (|< 7 i - <li\ < k < q ± + q 2 ) 

Jo 

x [2(N - 1) Pn{u)i,qi)pn(v 2 ,q 2 ) + 3 • 3! (uq, <71 )p CT (w 2 , g 2 )] , (28a) 

Imn(w, k; a) = Imll d (u,k;a) 

( A \rr \ ^ i i f 00 

—J 2(2n) 3 fc J dw i dw 2[l + /K) + /(w 2 )](5(a;-a;i-a; 2 ) 

POO 

x / dgi gi dg 2 g 2 0 (|<?i - g 2 < fc < <?i + 92 ) 

Jo 

x 4p cr (wi,gi)p 7I .(w 2 ,g 2 ) , (28b) 


where /(w) = [exp(w/T) — 1] _1 is the Bose-Einstein distribution function. These expressions are 
not ultraviolet divergent and thus do not need to be renormalized. 

For the real parts, we employ the following approximation. We only take into account the Hartree- 
Fock contributions, arising from the tadpole diagrams Figs. Hb, c (for the er-meson) and Figs. [3] a, 
c (for the pion). As discussed above, these contributions do not depend on energy and momentum. 
The real parts of the diagrams in Figs. [3 d, e and 0 d are functions of energy and momentum 
and are much harder to compute, involving an additional numerical integration as compared to the 
respective imaginary parts iP|). This has been done in the Hartree (or large-JV) approximation 
to the 0(N) linear sigma model in Ref. |.'17l |. In this approximation, the only energy-momentum 
dependent contribution arises from the diagram in Fig. 0 d. It was found in Ref. |3]| that this 
contribution leads to minor quantitative, but not qualitative, changes of the spectral density of the 
cr-meson. Therefore, in the present first step where we focus exclusively on the effects of a nonzero 
decay width, we neglect the real parts of the diagrams in Figs. [3d, e and [3d. Consequently [see 
Appendix 1X1 Eq. l1A81) |. 


ReS(cr) = ReE fc (n)+ReS c ((r) 

4A 2 r°° 

= Tr 7o""V3 / dudqq 2 [1 + 2 f(u)\ [3p CT (w, q) + (TV - 1) p w (w, g)] , (29a) 

A {2ny Jo 

Rell(cr) = Ren a (cr)+Ren c (cr) 

4A 2 r°° 

= ~N 0V)3 J dujdc l c l 2 [^ + 2f(u)} [p <7 (u>,q) + (N+l)p 7T {uJ,q)} . (29b) 


Since the real parts do not depend on energy and momentum, they only modify the masses of the 
cr-meson and the pion, cf. Eq. (33). The thermal parts of the self-energies in Eqs. (33b ~ /(w) are 
ultraviolet finite, but the vacuum parts, ~ 1, are ultraviolet divergent and have to be renormalized. 
The O(N) linear sigma model is an effective theory which is not valid above a certain energy or 
momentum scale. Therefore, we choose a particularly simple renormalization scheme, regulating the 
ultraviolet divergent part of the 10 integral in Eqs. (33 by introducing a cut-off A, cf. Appendix 
El In principle, we could also have chosen to regulate the q integral. However, the spectral density 
relates the energy to the momentum integral in Eqs. ( 33 . For large momenta q, the spectral density 
is of appreciable magnitude only in a range of energies lo ~ q, cf. Figs.[3and[3 Of course, our results 
will depend on the cut-off A; this dependence will be discussed in Sec. El The cut-off prescription 
is also used to compute the tadpole and the sunset contributions in Eq. ll(ial) for the condensate, 
cf. Appendix For a more formal treatment of renormalization within the CJT formalism, see for 
instance Refs. [l9l [30, [tH . 
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After obtaining real and imaginary parts HU and 03, respectively, one inserts them in the 
expressions m (if the imaginary part is zero) or 03 (if the imaginary part is nonzero) for the 
spectral densities. These spectral densities can then be used to again evaluate the real and imaginary 
parts of the self-energies. This defines an iterative scheme to self-consistently solve for the spectral 
densities as functions of energy and momentum. A convenient starting point for this scheme is the 
Hartree-Fock approximation, i.e., neglecting the imaginary parts altogether. 

The self-consistently computed spectral densities obey a sum rule j38|, 

/ OO 1 

^Ar,7r(w,k) = 1 . (30) 


In our calculation, there are two reasons why this sum rule may be violated. One is because we 
neglected the real parts of the self-energy diagrams in Figs. EJd, e and E| d. (The results of Ref. [37] 
suggest, though, that this has minor impact on the sum rule.) The other is due to the numerical 
realization of the above described iterative scheme. Numerically, one has to solve for the spectral 
density on a finite, discretized grid in energy-momentum space. If the imaginary part of the self¬ 
energy becomes very small, the spectral density converges towards a delta function. The support of 
the delta function may be located between grid sites, which causes a loss of spectral strength, which 
in turn violates the sum rule. 

Our prescription to restore the validity of the sum rule is the following. We first check whether 
the imaginary part is smaller than the grid spacing (in our calculations, 10 MeV in both energy and 
momentum direction) at the location of the quasiparticle mass shell, to = Lo a ^(k). If this is the case, 
we add a sufficiently wide numerical realization of the delta function, <5 num , to the original spectral 
density, p(oj, k) —► p(u>, k) + c ■ 5 num [u; — ui a ^(k)\, where c is a constant that is adjusted so that the 
sum rule is fulfilled on our energy-momentum grid, 



7— Wp CT ,7r(w, k) = 1 , 

Z7T 


(31) 


where <x> max is the maximum energy on the energy-momentum grid. 

On the other hand, if the imaginary part turns out to be sufficiently large, we presume that a 
possible violation of the sum rule <EB arises from neglecting the real parts of the self-energy diagrams 
in Figs. Eld, e and|3]d. In this case, we multiply the spectral density by a constant, p —► d ■ p, where 
d is adjusted so that the sum rule <EB is fulfilled. (By comparing the results to the case d = 1, we 
found that this somewhat ad hoc correction procedure does not lead to major quantitative changes, 
in agreement with the results of Ref. [371].') 

Restricting the sum rule to a finite range in energy as in Eq. m, however, causes the following 
problem. If the decay width of the particles is very large and consequently the spectral density a 
rather broad distribution around the quasiparticle mass shell, there will be parts which lie outside 
the energy-momentum grid. We could estimate the magnitude of this physical effect if we knew the 
behavior of the spectral density at energies to > w ma x. This is possible at zero temperature (with 
the help of Weinberg’s sum rules), but not at nonzero temperature. Here, we simply assume that 
this effect is sufficiently small to be neglected, i.e., we assume that a possible violation of the sum 
rule m is due to the two above mentioned artifacts and accordingly perform the correction of the 
spectral density. 

Another physical effect which causes a loss of spectral strength is if the quasiparticle energy 
Wo- iW (fc) > w max . In this case, we do not perform the correction of the spectral density as described 
above. This occurs at large momenta k, close to the edge of the energy-momentum grid. We checked 
that, in the numerical calculation of the integrals in Eqs. (EEJ, the integrands become sufficiently 
small in this region, so that the imaginary parts are not sensitive to this effect. 

This concludes the discussion of our improvement of the Hartree-Fock approximation by self- 
consistently including nonzero decay widths. 
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IV. RESULTS 


In this section we present the numerical results for the 0(4) linear sigma model obtained in the 
Hartree-Fock approximation improved by including nonzero decay widths, as discussed in Sec. urn 
The first step is to determine the parameters of the model. This is discussed in detail in Appendix 
m The zero-temperature values for the mass of the pion and the scalar condensate are M n = 
139.5 MeV and a = = 92.4 MeV ,36]. The cr-meson has a large decay width because of its 

decay into two pions. We take M a = 600 MeV and determine the width self-consistently within 
our improved Hartree-Fock approximation. For the cut-off required to renormalize the ultraviolet- 
divergent contributions of loop integrals, we consider two choices, A = 0 and 200 MeV. (As discussed 
in Appendix El it is not possible to consider arbitrarily large values for A; for M a = 600 MeV the 
maximum value is A max = 260 MeV.) For A = 0, we obtain for the parameters of the 0(4) model 
the values H = (121.60MeV) 3 , A = 19.943, and /r 2 = —(388.34MeV) 2 , while for A = 200 MeV, we 
have H = (107.16 MeV) 3 , A = 20.303, and /r 2 = -(411.38MeV) 2 . 

The condensate cr is shown in the left part of Fig. E as a function of temperature for the stan¬ 
dard and the improved Hartree-Fock approximation. The qualitative behavior is similar in the two 
approximations: the condensate drops significantly with temperature indicating the restoration of 
chiral symmetry. Since we consider the case of explicit 0(4) symmetry breaking by taking H ^ 0 
in Eq. 0 , the chiral phase transition is a cross-over transition. Nevertheless, one can define a 
transition temperature, T x , as the temperature where the chiral susceptibility da/dT assumes a 
maximum. Quantitatively, the inclusion of a nonzero decay width lowers T x as compared to the 
standard Hartree-Fock approximation by about 20% for A = 0, and by about 25% for A = 200 MeV, 
respectively. The transition temperature Tk — 160 — 175 MeV agrees well with recent lattice results, 
T x ~ 172 ± 5 MeV for the two-flavor case |4j. (Note, however, that the latter value is extracted from 
an extrapolation to the chiral limit, while our result is for the case of explicit symmetry breaking, 
i.e., for nonzero quark masses.) 

For the solution of the condensate equation ll(ial) the relative magnitude of the contribution 
from the sunset diagrams, Eq. (TeB, is negligibly small, of order ~ 10 4 , and thus can be safely 
neglected. In turn, not having to compute the integrals pertaining to the double loop, cf. Eq. EH, 

considerably speeds up the computation. _ 

In the right part of Fig. E we show the effective masses M a = m 2 + Re E and M n = 
\/ to 2 + Re n of the cr-meson and the pion as functions of temperature in both approximation 
schemes. Since the decay width of the pion remains comparatively small, cf. Fig. E the mass 
of the pion does not change appreciably when taking the nonzero decay width into account. On the 
other hand, the large decay width of the cr-meson at temperatures below T x , cf. Fig.E does influence 
the mass: in this range of temperatures the mass exhibits a stronger decrease with temperature in 
the improved Hartree-Fock approximation. At large temperatures T > T x , the decay width of the 
cr-meson becomes negligibly small, cf. Fig. E and the mass approaches the value computed in the 
standard Hartree-Fock approximation. Both cr-meson and pion masses approach each other above 
T x , indicating the restoration of chiral symmetry. 

In Fig. Ewe show the decay widths of cr-mesons and pions, defined as 13811391 


r a(k) 


Im T,[ui a (k), k ; cr] 

w ff (fc) 


r n (k) 


Im n[w 7 r (fc), k; cr] 

Un(k) 


(32) 


where u> a ^{k) is the energy on the quasi-particle mass shell, cf. Eq. (1221) . At small temperatures, 
due to the possible decay of a cr-meson into two pions, the decay width of the former is large, of the 
order of its mass. Note that the values of r CT obtained here at T = 0 are completely determined by 
the parameters to 2 , A, and H of the 0(4) model, i.e., without adjusting any additional parameter. 
They are reasonably close to the experimentally measured value r CT 600 - 1000 MeV 0. The 
cr decay width increases with temperature up to a maximum at a temperature T ~ 100 MeV and 
then decreases rapidly. The decay width of the pion vanishes at T = 0. It increases at nonzero 
temperature and assumes a maximum around T ~ 150 MeV. At this temperature the decay width is 
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about half of the corresponding pion mass. It decreases rapidly at higher temperature. Although the 
decay widths of both particles decrease at high temperatures, they do not become degenerate, the 
decay width of the cr-meson remains about a factor of 8 larger than that of the pion. This difference 
can be traced to the symmetry factors multiplying the one-loop self-energies for the cr-meson and 
pion, cf. Eqs. ( 1251 ) . For the a- meson, there is a factor 2 (TV — 1) in front of the pion loop and a factor 
3 • 3! in front of the cr-loop, so that the overall factor is ~ 2(TV — 1) + 3 ■ 3! = 24. For the pion, 
there is only a mixed er-pion loop with a symmetry factor 4. From this simple argument one already 
expects that the decay width of the cr-meson is about a factor of 6 larger than that of the pion. The 
remaining difference comes from the fact that the self-consistently computed spectral densities of 
tr-meson and pion under the integrals in Eqs. G9 are also different, cf. Fig. El 

One might argue that, at asymptotically large temperatures, effects from chiral symmetry breaking 
can no longer play a role and the decay widths, as well as the spectral densities, of cr-meson and 
pion should become degenerate. This is true in the chiral limit, where there is a thermodynamic 
phase transition between the phases of broken and restored chiral symmetry, and where cr = 0 in the 
latter phase. Here, however, we consider the case of explicitly broken chiral symmetry, where a > 0, 
even when the temperature is very large. Since the decay widths are proportional to cr 2 , they also 
do not vanish at large temperature. The difference in the symmetry factors for the self-energies of 
cr-mesons and pions then leads to different values for the decay widths and spectral densities. 

The self-consistently calculated spectral densities of the cr-meson and the pion as functions of 
the external energy u> and momentum k are shown in Fig. [5] and Fig. 0 for different temperatures, 
T = 80, 160, 240, and 320 MeV. Here we restrict ourselves to the case A = 0; the results for the case 
A = 200 MeV are quantitatively rather similar. For a detailed discussion, let us fix the momentum at 
k = 325 MeV and consider the spectral densities as functions of energy u> for different temperatures, 
cf. Fig. |5| 

At all temperatures, the pion spectral density (dotted line) exhibits a pronounced peak on the mass 
shell, at u>Tr(k) = yjk 2 + M% (cr). When the temperature is below 200 MeV, such that M n ~ 139.5 
MeV, cf. Fig. EJ the peak is located at w^(325MeV) ~ 350 MeV. Above T ~ 200 MeV, M n 
increases significantly with temperature, and the position of the peak shifts towards larger energies, 
04 . (325 MeV) ~ 500 MeV. The broadening of the peak is due to scattering of the pion off cr-mesons 
in the medium. 

In contrast to the pion spectral density, for temperatures be low ^170 M eV the a spectral density 
(full line) does not exhibit a peak at the mass shell, u> a (k) = i/fc 2 + Af 2 (cr). The reason is that u> a (k) 
is still sufficiently large to allow for the decay into two pions. Consequently, in this temperature range 
the cr spectral density is very broad. On the other hand, for temperatures above ~ 170 MeV, where 
oj a (k) drops below 2 M n , the two-pion decay channel is closed and the cr spectral density develops a 
distinct peak, whose width is due to scattering of the cr-mesons off pions and other cr-mesons in the 
medium. 

Two other features of the spectral densities shown in Fig.EJare noteworthy. The first is the region 
below the light-cone, K 2 = w 2 — k 2 < 0, where the mesons are Landau-damped. The second is the 
two-particle decay threshold. For cr — > 7T7t, this threshold is located at u> ~ 2 M n , for cr —► crcr, it 
is at u) ~ 2 M„ , and for 7r —■> cr7r at u> ~ M a + M T . The threshold is most clearly seen at large 
temperatures, e.g. T = 320 MeV, when both particles become degenerate in mass, M a ~ ~ 400 

MeV, and the threshold is at w ~ 900 MeV. 

V. CONCLUSIONS 

In this paper we have studied the O(N) linear sigma model at nonzero temperature within a 
self-consistent many-body resummation scheme. This scheme extends the standard Hartree-Fock 
approximation by including nonzero decay widths of the particles. In the Hartree-Fock approxima¬ 
tion, the self-energies of the particles consist of tadpole diagrams which have no imaginary part. 
Consequently, all particles are stable quasi-particles. In order to obtain a nonzero decay width, one 
has to include diagrams in the self-energy, which have a nonzero imaginary part corresponding to 
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decay and, in a medium, scattering processes. 

In order to incorporate the nonzero decay width in a self-consistent way, we apply the Cornwall- 
Jackiw-Tomboulis formalism. The standard Hartree-Fock approximation is obtained by considering 
only double-bubble diagrams in the 2PI effective action, leading to the (energy- and momentum- 
independent) tadpole contributions in the 1PI self-energies. In order to extend the Hartree-Fock 
approximation, we additionally take into account diagrams of sunset topology in the 2PI effective 
action. This has the consequence that the 1PI self-energies obtain additional energy- and momentum- 
dependent one-loop contributions which have a nonzero imaginary part. The spectral densities of 
tr-mesons and pions are then computed as solutions of a self-consistent set of Dyson-Schwinger 
equations for the er-meson and pion two-point functions, coupled to a fix-point equation for the 
chiral condensate. We only take into account the imaginary parts of the new one-loop contributions. 
We made sure that the spectral densities obey a standard sum rule by adjusting their normalization, 
if necessary. 

We found that the temperature T x for chiral symmetry restoration is about 20 — 25% smaller as 
compared to the Hartree-Fock approximation when including nonzero particle decay widths. Our 
value for the temperature of chiral symmetry restoration T x agrees reasonably well with lattice 
results 0. We computed the decay widths of cr-mesons and pions as a function of temperature. The 
vacuum value for the a decay width comes out to be in the experimentally observed range, without 
adjusting a parameter of the model. It stays approximately constant up to temperatures ~ T x and 
then decreases sharply with temperature. The pion decay width grows from zero at T = 0 to a value 
~ 100 MeV at T ~ T x , and then also decreases with temperature. 

We also investigated the spectral densities of cr-mesons and pions as functions of energy and 
momentum for temperatures in the range of T = 80 to 320 MeV. Below the chiral phase transition, 
the spectral density of the cr-meson is broad, due to the possible decay into two pions. It develops 
a peak at the quasi-particle mass shell above the chiral phase transition, when this decay channel 
is closed. On the other hand, the spectral density of the pion always exhibits a distinct peak at the 
quasi-particle mass shell. The width of this peak is due to scattering off cr-mesons in the medium. 
Above the chiral phase transition, the spectral densities of cr-mesons and pions become degenerate 
in shape. 

The present study can be continued along several lines. An important task is to include also the 
real parts of the self-energy contributions from the sunset diagrams. It would also be of interest to 
study the chiral limit, in order to see how the sunset diagrams affect the order of the chiral phase 
transition. In the Hartree-Fock approximation, the chiral phase transition turns out to be of first 
order, while universality arguments predict it to be of second order. 

Other possible future projects are the inclusion of other degrees of freedom, for instance strange 
scalar and pseudoscalar mesons, baryons jljj, and vector mesons The latter is of 

particular importance, since in-medium changes in the spectral properties of vector mesons are 
reflected in the dilepton spectrum |44| which, in turn, is experimentally observable in heavy-ion 
collisions at GSI-SIS, CERN-SPS and BNL-RHIC energies. 
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APPENDIX A: THE CALCULATION OF THE DIAGRAMS 

In this Appendix we discuss the calculation of the three types of diagrams (see Fig. 0 contributing 
to the self-energies and the condensate equation. To evaluate them explicitly, we use standard 
techniques of thermal field theory, see for example Refs. l38l llTl llfll. 
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1. The tadpole diagram 


In the following we discuss the calculation of the tadpole diagram in Fig. M a as a functional of 
the spectral density. In the imaginary-time formalism it is given by 

TST Y. (Al) 

n ' ' 

where T is the temperature, and A is either the a-meson or pion propagator. 

The first step is to perform the sum over the Matsubara frequencies. To this aim we introduce 
the mixed representation of the propagator 

,1/T 

A(-iu; n ,q) = / dr exp(-i^ n T)A(r, q), (A2) 

Jo 

with 

/ OO 7 

2 ^ [0(t) + f(u>)] p(u>, q) exp(-wr) , (A3) 

where fito) = [exp(w/T) — l] -1 is the Bose-Einstein distribution function. With the identity 


OO 

T exp (~iu! n T) = Y 


(A4) 


m— — oo 


the Matsubara sum can be performed analytically 

d 3 q rV T 


t = T Yf yy j 0 dTex p(- w « r ) A ( T -q) 

/' d 3 q 


ci /T 


(2tt) 3 


drA(r, q) 


5{t) + 6 r — 


T 


d 3 q 1 
2 

A(0,q) 


(27r) 3 2 
d 3 q 
(2tt) : 


A(0,q) + A ( -,q 


(A5) 


where use has been made of the KMS condition A(r, q) = A(r — 1/T, q). With Eq. IIA3I) we finally 
have 


T = 


dio f d 3 q 


2n J (27r) 3 


1 


/(^) 


p{u, q) ■ 


(A6) 


Due to isotropy of space, the spectral density of a scalar particle cannot depend on the direction of 
q, thus the angular integration can be carried out: 


^ , , dio r°° 

T = Att / —— 


dq 


2tt J o (27r) 3 


^ + /M 


P(v,q) ■ 


(A7) 


Using the fact that the spectral density for bosonic degrees of freedom is an odd function of the 
energy, p{u) = -p(-ui), 


T = 


4n l 

41 i ^(^ ?![i+2/m 


(A8) 
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This integral can be split into an ultraviolet-finite temperature-dependent part, 

TT = i ’l ^ (^ ?22/( “ ) ' ,( “’ ,) ’ <A9) 

and an ultraviolet-divergent vacuum part, which is renormalized by introducing a cut-off A for the 
lo integral, 

^ = 4?r L 3 q 2 9( ' A " ^ P ^’ ^ ’ ( A10 ) 

where T = T T + T A . 


2. The cut sunset diagram 


In this section we calculate the imaginary part of the cut sunset diagram of Fig. H b. In the 
imaginary-time formalism the cut sunset diagram is given by 

C(-iu>m, k) = [ Ai(/v - Q) A 2 (Q) = tV / ^-^Ai(-i(w m - w n ), k - q)A 2 (-iw„, q) (All) 
•/Q n J ( Zn ) 

where Ai,A 2 are the propagators of cr-mesons and/or pions. The major difference between the 
tadpole diagram, discussed in the last section, and this diagram is that it explicitly depends on the 
external four-momentum, = (—?'w m ,k). 

Analogously to the last section one introduces the mixed representation [see Eq. lA2l ] of the 
propagators Ai and A 2 

■t—v f d 3 q r 1 / T 

k) =T22 J ^ 277 j 3 J dnexp[-i(a; m -a; n )ri]Ai(Ti,k-q) 


rl/T 

/ dr 2 exp(-ia; n r2)A2(r2,q) . 
Jo 


(A12) 


The summation over the Matsubara frequencies leads to a delta function [cf. Eq. (IA4[l ] 

f d 3 q f 1 ^ T 

C(-iuj m , k) = J J drexp(-*w m r)Ai(r,k - q)A 2 (r, q) . (A13) 

Introduction of the spectral densities [see Eq. El] for both propagators in the mixed representation 
leads to 

C{-iuj m , k) = J ^ 3 J dr exp(-iui m T) J -^y[l + /(wi)]pi(wi, k - q) exp(-wir) 


diU2 
2i r 


[1 + /(w 2 )]/ 92 (w 2 ,q)exp(-w 2 r) . (A14) 


The integration over r can be performed analytically. Employing the identity 
c*h + w 2 ' 


exp - 


T 


- 1 


[1 + /( wi)][l + /( w 2 )] = -1 - /( wi ) - /( w 2 ) 


the expression can be rewritten as 

dw\ du -2 1 + /(wi) + /(w 2 ) f d 3 q 


r, ■ ^ P 

c(-zw m ,k)=y — 


27t -f- uq -(- tu 2 


(2tt) s 


pi(wi,k - q)p 2 (w2,q) ■ 


(A15) 


(A16) 
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Using the fact that the spectral densities do not depend on the direction of momentum and 

1 


d 3 x 2 t r 

f(\y-x\) 9 {x) = 


(2tt) 


V (27r) 3 


dx i x 1 dx 2 x 2 f(x 1 )g(x 2 ) 


x 0 (|xi - x 2 \ < y < Xi + x 2 ) , 


(A17) 


one can perform one angular integration, yielding 


C(-i(Jm, k) = 


X 

X 



duJi dm2 1 + f(wi) + f{m 2 ) 
27 t 27 t iui m + uj\ + ui 2 


2 tt 1 
k (27r) 3 
0 (\ qi - 


dqi qi dq 2 q 2 pi(uii, qi)p 2 (ui 2 , q 2 ) 


Jo 

921 < fc < 91 + 92) • 


(A18) 


The imaginary part of the retarded cut sunset diagram can be extracted using analytical continua¬ 
tion, —iuJm —> ui + ie, and the Dirac identity Iml/(:r + ie) = — 7r6(x ), 

1 1 f°° 

ImC(w,k) = p-- J duiidui 2 [l +f(uii) + f(ui 2 ))d(ui - uii~ ui 2 ) 

nOO 

x / dqiqidq 2 <?2©(|<7i -£2! < * < < 7i + <72) Pi(^i, <?i) £2(^2, <72) • (A19) 

Jo 

This integral is finite and therefore does not require renormalization. 


3. The sunset diagram 


In the imaginary-time formalism the sunset diagram shown in Fig. El c is given by the following 
expression: 

5 = T 2 ^ [ d~T^^.i(—iui m , 1)A 2 (—q)A 3 (— i(—— ui m ), — q — 1) . (A20) 

( 27r ) ( 2?r ) 

We introduce the mixed representation, see Eq. iE3>, for the three propagators Ai, A 2 , and A 3 . One 
can perform the Matsubara sums employing Eq. IA 41) . In the mixed representation one introduces 
the spectral densities pi,p 2 , and p 3 in order to perform the r integration analytically. Employing 
the identity 


exp 


U>1+U1 2 + Ul 3 

T 



[1 + /(wi)] [1 + f{ui 2 )} [1 + f(ui 3 )] 


= -1 - /(<*>l) - /(w 2 ) - f(w 3 ) - f(uii)f(ui 2 ) - /(wi)/(w 3 ) - f(ui 2 )f{u 3 ) , 


(A21) 


for the Bose-Einstein distribution functions we obtain 

3 duii d^2 du 3 1 + f{vi) + /(^ 2 ) + f(u 3) + /(m)/ 0 * 7 2 ) + f(vi)f{w 3 ) + f{ui 2 )f{ui 3) 
„ 27T 27T 27T 

d 3 Z d 3 ? 


5 = 


U7l + W 2 + ^3 

■pi(wi,l)p 2 (w 2 ,q)P3(w3 ,-q-l) . 


(A22) 


7 ( 27 r ) 3 ( 27 t ) 3 

Using the fact that the spectral densities do not depend on the direction of momentum, the angular 
integrals can be carried out with the result 

2 


5 = 


/ OO /*oo 

duii dui 2 dui 3 / dqi qi dq 2 q 2 dq 3 q 3 © (|9i - q 2 | <93 < 91 + 92) 
-00 J 0 


(2tt 

1 + f{uii) + f{ui 2 ) + /(w 3 ) + f(un)f{u) 2 ) + f(uii)f{ui 3 ) + f(ui 2 )f{ui 3 ) 


0>1 + W 2 + W3 

X Al(wi, 9 i)p 2 (w 2 ,9 2 )/93 (w 3,93) ■ 


(A23) 
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Using the antisymmetry of the spectral densities and relabelling integration variables, we finally 
arrive at 


5 = 


(2tt)' 


pOO 

■ / du i du >2 du > 3 dq i 91 dq 2 q 2 dq 3 q 3 0 (|gi - q 2 \ < q 3 < qi + q 2 ) 

Jo 

/(m)/(^ 2 ) + /(^l)/(^3) + /(^3)/(^2) + /(^l) + /<A 2 ) + /(^ 3 ) + 1 

<X >1 + W 2 + W 3 

/(wi)/(w 3 ) - f(ui)f(u> 2 ) + f(u> 2 )f(ui 3 ) + f(u> 3 ) 


UJ 3 U) 2 — U) 3 

/(^l)/(^2) - f{ui)f{u 3 ) + f(v2)f(v3) + /M 

Wl — W 2 + W3 

/(wi)/(w 2 ) + /(wi)/(w 3 ) - /(w 2 )/(w3) + /(wi) 


+ 

+ 

—+ W 2 + W 3 
X Pi(wi,gi)p 2 (w 2 ,g 2 )p 3 (w 3 ,g 3 ) • 

This integral can be split in an ultraviolet-finite temperature-dependent part, 
4 


(A24) 


= 


(2tt)' 


pOO 

■ / dwi (iw 2 dw 3 d<?i 91 dq 2 q 2 dq 3 q 3 0 (|gi - q 2 \ < q 3 < q 3 + q 2 ) 
Jo 


f{uJl)f( 0 J 2 ) + f(u>l)f(w 3 ) + f(u 3 )f{uj 2 ) 

OJl + U ) 2 + U ) 3 

/(^i)/(^ 3 ) - f{wi)f(u> 2 ) + f(u} 2 )f(u 3) 

Wl + W 2 — W3 

f(vi)f(u 2 ) - /(Al)/ 03 ) + f(u 2 )f(w 3 ) 

OJi - u 2 + w 3 

, /(wi)/(w 2 ) + /(wi)/(w 3 ) - f(uj 2 )f(u 3 ) 


—U)\ +U 2 +UJ 3 

X Pi(wi, qi)p 2 ( 0 J 2 , q 2 )p 3 {u> 3 , q 3 ) , 


(A25) 


an ultraviolet-divergent temperature-dependent part, which is renormalized by introducing a cut-off 
A for one of the 09 integrals, 


S A = 


(2tt) 7 


pOO 

/ dwi dw 2 dw 3 d<Zi <?i ^92 92 dq 3 q 3 0 (|<?i - 92 1 < 93 < 9i + 92 ) 
Jo 


/(wi)0(A - w 2 ) + /(w 2 )0(A - wi) + /(w 3 )0(A - wi) 


wi + w 2 + W 3 

/(w 3 )0( A-wi) /(w 2 )0(A-w!) /(wi)0(A-w 2 ) 


Wi + W 2 - W 3 Wi — W 2 + W3 

x pi (wi, 91 )p 2 (w 2 ,q 2 )p 3 (oj 3 , 93 ) , 


—Wi + w 2 + W3 


(A26) 


and an ultraviolet-divergent vacuum part, which is renormalized by introducing a cut-off A for two 
of the u>i integrals, 


<S AA = 


(2tt)' 


p OO 

■ / dwi dw 2 dw 3 dqi q 3 dq 2 q 2 dq 3 q 3 0 (I 91 - 9 2 | < 93 < 9i + Q2) 
Jo 


0(A — wi)0(A — u) 2 ) 
u>i + u> 2 + U) 3 


pi(wi, 9i)p 2 (w 2 ,92)P3(^3, 93 ) 


(A27) 


where S = S T + S A + S AA . 
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APPENDIX B: THE RENORMALIZED PARAMETERS 


In this appendix, we discuss how to determine the parameters /i 2 , A, and H of the O(N) model. 
At T = 0, the equations for the cr-meson and the pion mass are, cf. Eq. (1251) . 


Ml = m 2 + ^[3/ 2 + 3T a + (1V--1)T a ] 
Ml = ^ + ^[fl + r A + (N + 1)T a ] , 


(Bla) 

(Bib) 


where the zero-temperature contributions of the tadpole diagrams, 7^, are defined in Eq. ifTfoli . 
Solving these equations for A and fj, 2 one obtains 


A = A(A) = 


N Ml-Ml 
“8 / 2 + TA - 71 A ’ 

J 7T 1 11 (J TV 


M 2 = 




Equation 111 6 all for the zero-temperature condensate, <r = / w , can be written as 

f 8A fl 


H = f n { Ml- 


N 


[2(IV-l)5 AA (T + 3!5i A ]| 


(B2a) 

(B2b) 

(B2c) 


where 5 AA T is the zero-temperature contribution of the sunset diagram consisting of two pions and 
one cr-meson, and 5 AA CT is the corresponding one consisting of three cr-mesons, respectively. These 
quantities are defined in Eq. (TOTIl . Since these contributions are small, we neglect them. 

The term ~ H in the tree-level potential breaks the residual 0(N — 1) symmetry of the vacuum 
explicitly and favors a certain value for the scalar condensate cr = /, r . The sign of this value is 
determined by the sign of H. For H > 0, we must have a > 0. In order to preserve this relation, 
we read off from Eq. (IB2cl that A is restricted from above by _/VM 2 /(8/ 2 ). From Eq. llB2al) . this 
constraint leads to a restriction on the possible values of the cut-off A. For a given mass of the 
cr-meson, the maximum value for the cut-off, A max , decreases with increasing M a . For instance, for 
M a = 400, 600, 800 MeV, one obtains A max = 340, 260, 220 MeV, respectively. A similar behavior 
for the cr-meson mass as function of the cut-off has been discussed in Ref. 01 ■ 
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FIG. 1: The set of two-particle irreducible diagrams defining the improved Hartree-Fock approximation for 
the O(N) model, see Eq. ifil . The diagrams a,b, and c are the double-bubble diagrams, and d and e are the 
sunset diagrams. A full line denotes the full propagator for the a-meson and a dashed line the full propagator 
for the pion. The four-particle vertex ~ A is represented by a square and the three-particle vertex ~ A a by 
a triangle. 



FIG. 2: The self-energy of the cr-meson. The diagrams b and c are tadpole contributions generated by 
cutting an internal er-line in the double-bubble diagrams b and c in Fig. Q The diagrams d and e are 
one-loop contributions generated by cutting an internal cr-meson line in the sunset diagrams d and e of Fig. 
0 Lines and vertices as in Fig. 0 


a c d 



FIG. 3: The self-energy of the pion. The diagrams a and c are tadpole contributions generated by cutting an 
internal pion line in the double-bubble diagrams a and c in Fig.0 The diagram d is the one-loop contribution 
generated by cutting an internal pion line in the sunset diagram d in Fig.0 Lines and vertices as in Fig.0 
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FIG. 4: The values for the condensate (left panel) and the effective masses of the er-meson and the pion, 
(right panel) as functions of temperature. The values are calculated in the Hartree-Fock approximation 
(dashed and solid lines) and in the improved Hartree-Fock approximation (symbols) as discussed in Sec. lllll 



FIG. 5: The decay width of the cr-meson, r CT = ImE (ui)/u>, and the pion, F^ = ImII(w)/a!, in the improved 
Hartree-Fock approximation as a function of the temperature at momentum k = 325 MeV. 
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FIG. 6: The spectral density of the <r-meson as a function of energy uj and momentum k at temperatures 
80, 160, 240, and 320 MeV. 
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FIG. 7: The spectral density of the pion as a function of energy w and momentum k at temperatures 80, 
160, 240, and 320 MeV. 
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FIG. 8: The spectral density of the cr-meson and pion as a function of energy u> at temperatures 80, 160, 
240, and 320 MeV. The momentum is k = 325 MeV. 



FIG. 9: The general topology of the tadpole diagram a, the cut sunset diagram b, and the sunset diagram 
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